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Abstract — The mean key variable to control the burner 
performance for safe and efficient operations is the amount 
of excess air or oxygen which flow to react with fuel. In 
accordance with Typical Draft Profile in a Natural Draft 
furnace the numerical results of computing fluid dynamics 
model of the air flow from the environment to the burner by 
adjusting the area of air box prototype openings will give 
the necessary information to adjust the performance of 
burner as optimum. The model of air box prototype is useful 
tool in applying similarity theory on another burners 
designs which depends the air box tools to calibrate the 
quantity of combustion airflow to the burner. 

Keywords — excess air , air box design, computing fluid 
dynamics, natural draft furnaces. 

I. INTRODUCTION 

In many burners of natural furnaces today there are 
insufficient means to adjusting the ratio of fuel to air, the 
key variable to control the burner performance for safe and 
efficient operations is the amount of excess air or oxygen 
which flow to react with fuel. While some excess oxygen is 
necessary to have safe operation, too much oxygen wastes 
fuel, energy, and heat. The Air box is one component of 
burner which is used to control the amount of air excess by 
modifying its air registers. So the numerical results of 
computing fluid dynamics of the air flow from the 
environment to the burner by the air box will give the 
necessary information to adjust the performance of burner 
as optimum. 

In result; the process engineer will not depends on the 
laboratory analysis about the flue gas data from all sides of 
the combustion chamber of furnace at each burner. That is 
because the numerical results give a lots of information 
about the air flow through the air box to burner and its 
quantities from a design point to the operating case as a 
result of coupling the average velocity of combustion air 
which enter from the gats of air box into burner with the 
pressure arise across the burner or the air box. 


The pressure arise across combustion equipment like a 
burner etc., can be calculated from the following equation 


[ 1 ] 



Pa 


U = velocity of air at burner exit m/s 
K is the flow resistance factor where K=1.5 for tangential 
burner and K= 3 for swirl burner of gas and heavy oil [1] 
But from operation point, the air flow rate in burner can 
calculated from [1] 


273 + T air 
V . = B V°(n ■ ) -— 

v cur u ' v \ ,L airJ' 273 

V air The flow of air m 3 /s 
V° Theoretical air flow Nm 3 /kg 
B Fuel feed rate kg/s 

n air Excess air coefficient at furnace outlet (there is not 
leakage into furnace) 

T air Air temperature °C 

So the numerical results information which are gotten from 
the prototype air box can be used to predict the quantity that 
may be enter from the doors of the air box by using the 
theory of dynamics similarity which related by the 
Reynolds number of the air flow into the burner at all times 
with its Reynolds number at design case. 

From a design point these gates must provide enough air for 
fuel’s burning process. The correlative with maximum 
amount fuel consumption B max and surface area of all gates 
F g is defined by [2] 


F g (m 2 ) 


Bmax ^ C 

Wg x / X K 


B ma x Maximum amount fuel consumption 
l a Amount of air providing for burning the fuel 
j Area of activated gate = 1 in this case study 
k Discharge consumption coefficient 0.6-0.75 
w g Air speed through gate 0.5 m/s 

Here the process engineer thinks that in fact the excess air 
does not related only by adjusting the air box registers 
because it is important to maintain the correct furnace draft. 
So too little draft can damage the metal structure and snuff 
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out the burners or flame and too much draft can pull 
excessive amounts of air into the furnace. So the operation 
of fired furnaces is complex process and there are different 
variables on both the combustion and process sides of the 
furnace. In result; there is just enough air to assure complete 
combustion. Once, the air registers are adjusted for proper 
combustion, the damper should be adjusted for proper draft. 
So the engineer can make natural draft Heater adjustment 
according to the following flow chart showed in fig (1). 

TO benefit from the numerical results information which 
are gotten from the prototype air box from a design point, 
the process engineer must begin to ultimate the pressure 
arise across the burner (draft available at burner level),and 
he must follow these general rule which are mentioned for 
any duct in furnace not only for burner . 



Fig. 1: Natural Draft Heater Adjustment Flow Chart [3] 


1- General Rules to ultimate the pressure arise across 
any air duct in furnace 

By adding the natural draft components of all parts of the 
air ducts, the total natural draft YP n d is obtained, so the 
total pressure drop A P air in the air ducts in the furnace can 
be obtained by subtraction [1] 

kPair = kP r — Y Pnd P a 

A P air Is the total pressure drop in the air ducts in the 
furnace, Fig.(2) Shows typical draft Profile in a natural draft 
furnace 

A P r Is the total pressure drop after correction of elevation. 


AP r = 


i 


101325 

-p— 

r av 


The sum of resistance of all parts of the air duct gives the 
total resistance of the air duct£ A P. 



pressure pressure 

(SE) r + R 0 ~Pa * - -* 

(SE) C + (SE) s = P b + P c + Stack exist toss ♦ R 0 
A P, * Pressure drop across burners (draft available at burner level) 

SP b = Flue gas pressure drop across convection section 
A P c = Losses due to damper, stack friction, stack entrance 

Fig.2: Typical Draft Profile in a Natural Draft furnace [3] 


H Pnd Is the total natural draft of all parts of the air duct 


- Natural draft in air ducts in furnace can be calculated from 
Pnd = H.g . (pa - Pa) Pa 

Pnd Is the natural draft in air duct. 

H Is the height between the start and end sections of the 
duct m 

g Is the Gravitational constant 9.8kg/s 2 

p a Is the density of air at the lowest point of duct kg/m 3 

p a Is the density of air at the highest point of duct kg/m 3 

2- Mathematical model on air flow into swirl burner as 
case study 

The burner used in this study consists of two concentric 
tube with a tangential and radial swirler placed in annular 
arrangement. Fig. (4) and (5, B) illustrates the important 
aerodynamics feature of the vertical low NOx swirl burner. 
The combustion air is introduced through two concentric 
annular nozzles of which the positions of the primary, 
secondary air, and fuel tube are adjusted to produce the 
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desired internally staged flame, the burner utilizes swirl 
induces centrifugal body forces to dampen turbulent 
exchange between the fuel and air stream allowing an 
extended residence for fuel conversion chemistry in locally 
fuel -rich environment prior to burn out [4] .The NOx 


emission is at its peak at stoichiometric conditions when 
other pollutants, like CO and VOC, is at its minimum (see 
Fig.(3)).This is mainly because of the strong relation 
between temperature and thermal NOx production [5] 
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Fig. 3: NOx formation rate driven by temperature at equivalence ratio A and Schematic drawing of (a) air staging , and (b) fuel 

staging [5] 


NOx reduction needs to apply correct operating parameters 
and used as a design basis, so in accordance with NOx 
reduction work, it is essential that excess air levels be 
reduced to close to 2-4% O 2 as all 

NOx emissions are corrected to3% O 2 (see fig. (4))[3]. The 
most significant useful phenomenon of the swirling flows is 
the existence of internal recirculation zone IRZ which aids 



in flame stabilization by providing a storage of hot flow flue 
gas, so its flow field is thank for its effect in internal flue 
gas recirculated to the combustion zone to reduce the NOx 
emissions, by reducing the oxygen concentration as well as 
lowering the flame temperature. In internal FGR, the fuel 
gas is recycled into the flame zone due to burner 
aerodynamics [5] 



Fig. 4: Excess Air vs. Oxygen Content in Flue Gas (A) [3] and typical low NOX swirl burners flow field with internal staging 

schematic (B) [4] 
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Fig. 5: Shows a schematic diagram of the air and fuel flow 

In result, combustion product are recirculated and by this 
phenomena reduced velocity region is created, the boundary 
of the recirculation eddy is determined by the radial points 
at which the forward mass flow equals the reverse mass 
flow at the axial station so the region which is enclosed 
where the axial velocity reaches zero magnitude in each 
axial plane [6]. 

3- Computing fluid dynamics for the air box of burner: 
-air box geometry: 

Air box Consists of Two cylinders, the internal one is fixed 
with thickness 2 mm and dimensions 804 1 x370 mm and 
diameter 258 ° 5 mm, and the outside is movable one with 
thickness 2 mm, dimensions 823 +2 x268 mm and diameter 
264 +1 mm. These two cylinders have the same gates and by 
moving the outer cylinder against the fixed cylinder, 
changes the air flow area. See (fig. (6) Air box fluid 
domain) 


velocities of a swirled jet burner (a) [6] and burner geometry (b). 



Fig.6: Air box fluid domain (full opining) 

Ambient air temperature, pressure and humidity affect the 
air flow rate if the air temperature increases the air flow rate 
decreases. As definite temperature, pressure and humidity 
of air increase, water vapor in ambient air is decreased and 
leads to decrease of the percentage of O 2 . So the air flow 
rate through the air box at definite pressure drop by using 
corresponding - state is defined by [7] 

Pair^b^Vb 

mair = J k p 


www.awublications.com/iicmp 


Page | 7 



















International journal of Chemistry, Mathematics and Physics (IJCMP) [Vol-2, Issue-1, Jan-Feb, 2018] 

https://dx.doi.Org/10.22161/iicmp.2.l.2 ISSN: 2456-866X 


m air Air mass flow rate kg/s, g Gravitational constant 
9.8kg/s 2 ,i4 d Cross section area in internal cylinder m 2 , Ap b 
Air side pressure drop across air box in pascal (natural draft 
is available at air box level),/c p Pressure loss coefficient 
through air box and can be estimated as equal to 0.6-0.75 
- Cold mathematical model flow of vertical air box: 

As in the case of any problem in flow filed this particular 
type of flow described by the mathematical expressions of 
continuity and momentum in Cartesian coordinates: 
Continuity equations: 

du dv dw 


dx dy dz 


0 


Momentum Equations: 

In this situation the body force with constant density equals 
to the balance between the driving force of buoyancy and 
the initial hydrostatic pressure, so at any point in the fluid 
domain the equation of motion is: 


du t d{ujU J 
dr dXi 


dp d T u 


Where: T is the actual temperature applicable locally. : 
is the undisturbed constant temperature. B T \ is the 
coefficient of the thermal expansion of air and by using the 
equation of stat for the perfect gas 

p 


P = 


RT 


Where R is universal gas constant, p is the local static 
pressure, and p is the local density so ( B T = 1/T ), T \j is 
the Reynolds stress tensor divided by the density and it is 
computed by using the Boussinesq model equation 


T ij = u'jiij = v T 


du t duj 
dxj dx t 


2du k \ 2 

3 ^]- 3 K ^ 


Where K is the turbulent kinetic energy, 8 t j is Kronecker 

Delta Function, Then the letters i, j, k ...can be used as 
variables, running from 1 to 3 so Instead of using x, y, and z 
to label the components of a vector, we use 1, 2, 3 and v x is 
the turbulent eddy viscosity can be calculated by K-Epsilon 
model so 



Where is constant equals to 0.09 
These models are not valid near wall to model wall effect so 
wall function have to be employed. 


- Wall y+ Strategy for Dealing with Wall-bounded 
Turbulent Flows for air box geometry: 

A wall y+ in the range of 30 to 60 is determined to be 
sufficiently accurate. Where the log-law region is resolved. 


It is also advisable to avoid having the wall-adjacent mesh 
in the buffer region since neither wall functions nor near 
wall modelling approach accounts for it accurately, so for 
the existence of deferent boundary layers in air box 
geometry enhanced wall function is used 

- Boundary condition: 

Here the air box was opened to the environment from 
surrounding sides of the gates at the maximum fuel gas 
consumption which it was B max = 0.02 Nm 3 /s and 

( l a = { ^^ air ) ) • I n result from design point, air mass 

flow will be equal to (p a ir- B max . l a ) at X = 1.15 
So calculating the air mass flow yields 

- 3 ———) ), and the following are the technical 

1Nm ful,X= 1 . 15 / 

data of burner: 

- Burner axis: vertical 

- Regulating range of capacity of burner 20%-60% 

- Excess of air nair 2=1.15 

- Required type of the burner: gas type with natural draft of 
air fitted with a gas injector stabilization torch; Hydrogen 
content in the heating gas must not exceed 60% volume, 
and the fuel doesn’t contain any nitrogen compounds 

- Chemical composition %-vol for gaseous fuel: 

H 2 -28, Ci- 34, C 2 -17, C 3 -17, C 4 -3, C 5 - 1, H 2 0-1.5, H 2 S- 
max.5mg/Nm 3 

* Net calorific value: llOOOkcal/kg 

* Density kg/Nm 3 0.924 

* Flammability limits % 3.4-16.5 

- Temperature in the spot of the burner inlet 30°C 

- Heating gas consumption of I burner Nm 3 /h; 

For min. output 23.2, for nominal output 58.0 and for max. 
Output 69.6 

II. RESULTS AND CONCLUSIONS 

The wall y+ is a non-dimensional number similar to local 
Reynolds number, determining whether the influences in 
the wall-adjacent cells are laminar or turbulent, hence 
indicating the part of the turbulent boundary layer that they 
resolve [8], and it can be seen below as well as Velocity 
vector plots, can be seen below in the photos of Fig. (a), (b) 
and (c). These plots give a lot of information about the flow 
field of the air box of the burner. The recirculation at the 
corners of air gates and air flow by the area of gates is also 
obvious from these figures. The air combustion tends to 
move into the air box by natural draft of the furnace at its 
level. 
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first parte 1 fem4siml 
Load Casa 1, Static 
Velocity - 
Mn : 0.000, 

I bftAM 

8.185 
, 7.275 
6.366 


5.457 
4 547 
3.638 

( 2.728 
1.819 
0.909 
0/000 

Units = m/sec 


first parte 1_fem4_*im1 : Solution 1 Result 
Load Case 1. Static Step 1 
Y Plus • Element-Nodal* UnatMted, Scalar 

f 7.40E+002 

271 65 m 

99.78 
36.65 
| 13.46 
4 95 

1.82 sC: 

■ 0.67 / 

■ 0 .^45 
"~*O90 
Units = Unitlass 



first parte t Jem4_sim1 : Solution 1 Result 
Load Case 1. Static Step 1 



Velocity - Element-Nodal, 
Min : 0.01, Max 

1 9094 28 
8184.85 
■ 7275.42 
6366 00 
5456.57 
4547 14 
3637.71 

1 2728 29 
1818 86 
9()9.43 
0 * 1 - 

Units ■ mm/sec 


first parte 1_fem4_sim1 : Solution 1 Result 


Load Case 1, Stabc Step 1^| 

Y Plus - Element-I^ottfc Unaveraged Si 

m. .■ 


I 7.40E+0I 
271 65 
99.78 
36.65 
13.46 
4 95 
1.82 
I 0.67 
■ *245 
■ 0*90 
Units = Unitlass 


i\ 1 


-M 


\ 




first parte 1_rem4_sim1 Solution 1 ResJt 
Load Case 1 , Static Step I 
Y Plus - Element-Nodal, 

Min . 0.09. Ma« : 2013.39. Units 

1 2013.39 
1812.06 
n 1610 73 
1409.40 
1208 07 
1006.74 
805.41 
604.08 
k 402.75 
I 201.42 
■~0 % 09 

Units = Unitless 


Fig. a: y+ and velocity vector plots in the field of air box for solution 1. 


Solution 1 by iteration method at a steady state-physical relaxation step equals to 0.05 sec, the convergence criteria of RMS 
residuals is 0.0002, first order is used for advection scheme, distorted fluid elements are 6927elements , and the following table 1 
ultimate the data of this solution 


Table. 1: Ultimates the data of solutionl 


Computing Reynold 
number 

Number of 

element 

y+ information on adiabatic faces 

Computing cost 

1.8xl0 5 

3810849 

min y+ 

max y+ 

ave. y + 

52min 24sec 

1.59x10-' 

3.28x10' 

5.72 


first parte 1 (1)_sim1 : Solution 1 Result 

Load Case 1, Static Step 1 

Velocity - Element-Nodal, Averaged, Magnitude 



Units = m/sec 


first parte 1 (1)_sim1 : Solution 1 Result 
Load Case 1, Static Step 1 
Velocity - Element-Nodal, Averaged, Magnitude 
Min ^0*00, Max : 9.152, Units = m/sec 

• 8.236 
7.321 
6.406 
5.491 
4.576 
3.661 
2.745 
u h830 
■ 0^915 
1 0.000 
Units = m/sec 




Units = Unitless 



first parte 1 (1)_sim1 : Solution 1 Result 
Load Case 1, Static Step 1 
Y Plus - 


284.14 


y 0.^95 
1 0.T55 


first parte 1 (1)_sim1 : Solution 1 Result 


Load Case 1, Static Step 1 



Units = Unitless 


Units = Unitless 


first parte 1 (1)_sim1 : Solution 1 Result 



Load Case 1, Static 
Y Plus - Elem 


7 

284.14 
111.04 
43.39 
16.96 
6.63 
2.59 
, 1.01 
| 0.395 
■jO.155 


Units - Unitless 


Fig.b: y + and velocity vector plots in the field of air box for solution 2. 
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Solution 2 by iteration method at a steady state-physical relaxation step equals to 0.05 sec, the convergence criteria of RMS 
residuals is 0.0002, first order is used for advection scheme, distorted fluid elements are 22432 elements , and the following table 
2 ultimate the data of this solution 


Table. 2: Ultimates the data of solution! 


Computing Reynold 
number 

Number of 

element 

y+ information on adiabatic faces 

Computing cost 

1.8xl0 5 

1236619 

min y+ 

max y+ 

ave. y + 

12min 20sec 

3.70x10'' 

4.08x10' 

6.26 


first parte 1 (1)_sim1 : Solution 2 Result 
Load Case 1, Static Step 1 
Velocity • Element-Nodal, Availed. Magnitude 
Min : 0800, Max 9.154, Units “ m/sec 

1 8.239 £ 

7.323 
6.408 
5.493 
I 4.577 
3.662 
2.746 
1.831 
0.915 


first parte 1 (1)_sim1 : Solution 2 Result 
Load Case 1, Static Step 1 
Velocity ■ Element-Nodal, Averaged. Magnitude 
Min : 0 000, Max 9 154, Units ■ m/sec ^ 


lH- 




Units = m/sec 

first parte 1 (1) siml : Solution 2 Result 
Load Case 1 Static 

Y Plus • Element-Nodal J^eraged Scalar 

Wf 003. un,ts. unit,... 

" 6.74E+00# ^ 

248.33 
91.48 
33.70 

! 12.41 r ~ 

4.57 

1 1.68 
0.62 

Qjm 

Units = Unitless 



first parte 1 (1)_sim1 : Solution 2 Result 
Load Case 1. Static Step 1 
Y Plus - Element-Nodal, Aven 
1 83E+003, l 




first parte 1 (1) siml : Solution 2 Result 
Load Case 1 , statitrajjb,, 

Y Plus - Element-Nodal, Averaged, Scalar’ 
. 1.83E+003, Units 



' 0.1384 

Units = Unitless 



first parte 1 (1) siml : Solution 2 Result 
Load Case 1, Static Step 1 . 

Y Plus - Element-Nodal, Averagod^Scalar 

Unit’s -Unitless 



Fig.c: y+ and velocity vector plots in the field of air box for solution 3. 


Solution 3 by iteration method at a steady state-physical relaxation step equals to 0.05 sec, the convergence criteria of RMS 
residuals is 0.0002, first order is used for advection scheme, distorted fluid elements are 331 elements , and the following table 3 
ultimate the data of this solution 


Table. 3: Ultimates the data of solution3 


Computing Reynold 
number 

Number of 

element 

y+ information on adiabatic faces 

Computing cost 

1.8xl0 5 

1922514 

min y+ 

max y+ 

ave. y + 

23min59sec 

2.78x10-' 

3.24x10' 

5..44 


The comparison between these three solutions emphasis that solution 2 is preferred. That is because it gives the necessary 
information by the best computing cost as it is clear from fig. (c) Which shows the same velocity profiles for the three solutions 
at the same path. 
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Fig.d: Shows the same velocity profiles for the three solutions at the same path. 



Fig.e: Which shows low NOx burner with two air boxes 


In fact, by industry these kinds of air boxes are used as tools 
to adjust the fuel /air ratio in many burners’ designs which 
work on natural draft as it can be seen from fig. (e) Which 
shows low NOx burner with physical separate between two 
stages (two air stages and two fuel stages). Here burner 
depends on two air boxes for natural draft of the 
combustion air. So in result this model for air box in the 
low NOx swirl burner which is studied is too useful to get 
the important information and calibrate another burners 
performances by using the similarity theory as dynamics 
and geometry parameters, or to build different air boxes 
models by the same methodology to get the best from the 
system as large efficiency and less pollution especially for 
adjusting NOX emissions. 
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